Vortex sheet in rotating two-component Bose-Einstein condensates 



o 
o 

(N 
D 



I 

O 

o 



> 

o 

(N 
00 

O 

X 
J3 



Kenichi Kasamatsu^ and Makoto Tsubota^ 
^Department of Physics, Kinki University, Higashi- Osaka, 577-8502, Japan 
'^Department of Physics, Osaka City University, Sumiyoshi-Ku, Osaka 558-8585, Japan 

(Dated: February 16, 2009) 

We investigate vortex states of immiscible two-component Bose-Einstein condensates under rota- 
tion througli numerical simulations of the coupled Gross-Pitaevskii equations. For strong intercom- 
ponent repulsion, the two components undergo phase separation to form several density domains of 
the same component. In the presence of the rotation, the nucleated vortices are aligned between the 
domains to make up winding chains of singly quantized vortices, a vortex sheet, instead of periodic 
vortex lattices. The vortices of one component are located at the region of the density domains 
of the other component, which results in the serpentine domain structure. The sheet configuration 
is stable as long as the imbalance of the intracomponent parameter is small. We employ a planar 
sheet model to estimate the distance between neighboring sheets, determined by the competition 
between the surface tension of the domain wall and the kinetic energy of the superflow via quantized 
vortices. By comparing the several length scales in this system, the phase diagram of the vortex 
state is obtained. 

PACS numbers: 03.75.Lm, 03.75.Mn, 05.30.Jp 



I. INTRODUCTION 

Diverse interests have been focused on the quan- 
tized vortices in trapped atomic Bose-Eistein conden- 
sates (BECs) flj. This pristine Bose-condensed system 
is highly controllable, having a great advantage to inves- 
tigate the physics of topological defects such as vortices. 
Rotation of the trapping potential can create a lattice of 
quantized vortices [3| , whose properties are observable by 
direct imaging. The next step of this kind of study will be 
elucidation of the vortex states in multicomponent BECs 
0. A multicomponent order parameter allows the exis- 
tence of the various unconventional topological defects, 
which have been discussed in the other condensed matter 
systems and cosmological theory 0, Q . 

Two-component BECs, the simplest example of the 
multicomponent condensates, have also attracted much 
interest to study the novel phenomena not found in the 
single component system. The typical example is phase 
separation of the binary BEC mixtures, where each com- 
ponent separates spatially due to the strong intercompo- 
nent repulsion. Several ex per iments have observed this 
phase separation [6l.[7l.[8l.[9l. ll0l . [Tl| , su pported by the nu- 
merous theoretical studies [11 111 El, [11, [l3. Ill III ■ 
We can now control miscibility or immiscibility of the 
two-component BEC, because intercomponent interac- 
tion is tunable by changing the s-wave scattering length 
via a Feshbach resonance, as demonstrated recently 

MM. 



Rotating two-component BECs could have a rich vari- 
ety of vortex phases in addition to the conventional tri- 
angular lattice i2d,[lll[23,[2l,[2l[2|]. Schwcikhard et al. 
observed the transient structural transition of the vortex 
lattice in two-component BECs [26'] . The richness of the 
vortex structure originates from the intercomponent in- 
teraction. For instance, with increasing the strength of 
intercomponent interaction from zero, the vortex lattice 



deforms from triangular to square, eventually evolving 
multiple vortex sheets with each coniponent made up of 
chains of singly quantized vortices (20l . [2l|. Therefore 
it is challenging to see experimentally these structural 
transitions by using two-component BECs with tunable 
interactions. 

In this paper, we study the detail of the vortex sheet 
structure in rotating immiscible two-component BECs, 
which has not been explored so much. This system has 
three interatomic coupling constants denoted by gi and 
(72 (for intracomponents), and 1712 (for intercomponent). 
We confine ourselves to the phase-separated regime; in a 
homogeneous system, the condition is given by [Tsl . [Tg} 



512 > 51.92- 



(1) 



To give a clear picture of the vortex state, we first address 
the simulation results of the coupled Gross-Pitaevskii 
(GP) equations, discussing various configurations of the 
vortex states in immiscible BECs. Next, we give an an- 
alytical discussion of the sheet structure. The distance 
of the vortex sheet can be estimated by the ratio of the 
surface tension of the domain boundary and the kinetic 
energy of the superflow due to the vortices [13] ■ Compar- 
ing the several length scales characteristic of our system, 
we can obtain the phase diagram of the vortex states, 
which provides observable conditions of the several vor- 
tex phases, such as multiple vortex sheets with stripe pat- 
tern, or "serpentine" sheets, or the "rotating droplets." 

There are other condensed-matter systems character- 
ized by the multicomponent order parameters in which a 
vortex sheet is observable [13, [II]- In Sec. |TT1 we give a 
brief review of a vortex sheet in supcrfluid "^He to help the 
following discussion. Section [1111 consists of three parts; 
after formulating our problem in Sec. IIII Al we present 
the numerical results in Sec. IIII Bl and analytical discus- 
sion in Sec. IIII CI for the prolem of two-component BECs. 
We conclude this paper in Sec. IIVI 
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FIG. 1: (Color online) Vortex sheet array suggested by Lan- 
dau and Lifshitz. Straight vortex lines are aligned along cylin- 
drical planar soliton. The schematic view of the radial velocity 
profile is also shown. 



distance b is determined by the competition of the sur- 
face tension a of the sohton and the kinetic energy of 
the counterflow Vn — Vg outside the sheet, where v„ is 
the normal fluid velocity. The planar sheet is assumed to 
be parallel to the x axis for a large system. The motion 
of the normal component makes a rigid body rotation 
V X Vn = 2f7z and can be represented by the shear flow 
v„ = — 2riyx parallel to the planes. In the region be- 
tween nearest planes, the vortex-free superfluid velocity 
Vg is constant and equals the average Vn to minimize 
the counterflow. Thus the velocity jump across the sheet 
is = 2flb. The counterflow energy per volume is 
(1/6) / ^Ps(vn - ^s)^dy = where ps = mus 

is the superfluid mass density. The surface energy per 
volume equals a/b. Minimizing the sum of two contri- 
butions with respect to b, one obtains b = (Scr/psil^)^/'^. 
For Q = 1 rad/s, this gives 0.3-0.4 mm, which was quanti- 
tatively conflrmed by the NMR absorption measurement 
in superfluid ■^He-A The possible sheet configura- 

tions and their response for rotating drive were studied 
by numerical simulations in Ref. |30|. 



III. VORTEX SHEET IN TRAPPED 
TWO-COMPONENT BECS 



II. VORTEX SHEET IN A ROTATING 
SUPERFLUID HELIUM 



The vortex sheet is well known from classical tur- 
bulence as a thin interface across which the tangential 
component of the flow velocity is discontinuous. In the 
limit of vanishing width of the sheet, the vorticity ap- 
proaches infinity within the sheet. Historically in su- 
perfiuids, an array of vortex sheets was suggested as the 
possible ground state of rotating superfluid *He [2^ . The 
schematic picture of the vortex sheet, suggested by Lan- 
dau and Lifshitz, is shown in Fig. [T] A vortex sheet 
consists of a chain of quantized vortices with the same 
orientation of circulation, which are confined by the cylin- 
drical planar soliton. The tangential component of the 
superfluid velocity Vg discontinuously jumps across the 
sheet. It soon turned out, however, that in superfluid 
^He a vortex sheet is unstable towards breaking up into 
separated vortex lines, namely forming a vortex lattice. 

In contrast, in superfluid "^He-A, the vortex sheet is 
stable due to the confinement of the vorticity within the 
topologically stable solitons [131 . In this system, vor- 
tices are confined by the planar soliton sheets, which 
separate regions with opposite orientations of the orbital 
angular momentum (1 vectors) of the Cooper pairs. In 
equilibrium, the vortex sheets uniformly fill the rotating 
container, forming accurately equidistant layers with the 
distance b. Here, the calculation by Landau and Lifshitz 
is useful to explain the equilibrium distance between the 
planes of vortex sheet in ^He-A in a container rotating 
with riz, even though their calculation did not take into 
account the quantization condition |,29i] . The equilibrium 



A. Coupled Gross-Pitaevskii equations 

We consider atomic two-component BECs, character- 
ized by the condensate wave functions 4"^ {i = 1, 2). The 
condensates are assumed to be confined by trapping po- 
tentials Vi{r) = mi{ui'^r^ -|-wfjZ^)/2; we neglect the shift 
of the potential minima for two components due to the 
gravitational effect or the difference of the hyperfine sub- 
levels. The potential is assumed to rotate at a rotation 
frequency about the z axis as $7 = flz. Viewed from 
the frame of reference corotating with the trap potential, 
the GP energy functional of our problem reads 



^ \-A — 1 o 



2y72 



2m i 



(2) 



Here, rrn is a mass of the ith atom, = iMl{ydx — xdy) 
is a centrifugal term, and the coefficients gi {i — 1,2) and 
9i2 represent the atom-atom interactions. These inter- 
actions are expressed in terms of the s-wave scattering 
lengths Oi and a2 between atoms in the same component 
and ai2 between atoms in different component as 



9i 



m,. 



912 



2TTh ai2 

mi2 



(3) 



where 



m^^ +m2^ is the reduced mass. Each com- 
ponent interacts with the other through the intercom- 
ponent mean-field coupling (?i2|^'ip|^2p, which yields 
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structures and dynamics not found in a single- component 

BEG @, Mill [11, [II, [13, [II, [11, IS E El- 

The properties of two-component BEGs are described 
by the coupled GP equations, obtained from Eq. ^ by 
using a variational procedure ihd'i/i = SE/d'i>* as 



ih 



dt 



2y72 



2to,- 



V^+g^\^,\^'+gl2\■^3^,\^ 



(4) 



Since we are interested in the stationary solutions of this 
equation, we consider the time-independent coupled GP 
equations by substituting 5'i(r, i) ~ \I>i(r)e~*'^''/^ as 



2m, 



(5) 



where we have introduced the Lagrange multiplier /i; 
which represents the chemical potential and is deter- 
mined so as to satisfy the conservation of particle number 

Even for = the equilibrium solutions of the cou- 
pled GP equations ([5]) exhibit a rich variety of struc- 
tures, depending on the various parameters of the system 
[I1,[I1[IJ,Ii1[i1,[I3,[I1- in particular, the intercompo- 
nent interaction gi2 plays an important role in determin- 
ing the ground state structure; when the intercomponent 
interaction is strongly repulsive, the two components 
phase separate. Adapting the Thomas-Fermi approxi- 
mation, one can see that there is no solution with over- 
lapping density profile when > 5152 [H, [H, [13, [Hi- 
This regime will be focused on in the following discussion. 
Then, in the space in which only one density is nonvan- 
ishing {rii = l^'ip 7^ and nj = 0, for i ^ j), the density 
profile is given by the single-component Thomas-Fermi 
profile Hi = [fii — Vi(r)]/gi. Since the Thomas-Fermi 
approximation fails to describe the domain boundary re- 
gion, more careful analysis is necessary to determine the 
explicit density profile in this region (sil . [s^ . 



B. Numerical results 

There are a few papers discussing the properties 
of immiscible two-component BEGs under rotation. 
Kasamatsu et al. revealed numerically that the wind- 
ing vortex sheets appears as a stable vortex state instead 
of the periodic vortex lattice 



2l|. Woo et al. studied 



Tkachenko excitations in rotating two-component BEGs, 
finding that some highly excited Tkachenko modes can 
give rise to the shear flow of vortices along the vortex 
sheet [2^. In a slowly rotating limit, Malomed et al. 
studied the stability of the rotating cross pattern of do- 
main walls separating one component in the first and 



third quadrant in the x-y plane from the other in the 
second and forth quadrant [33| . 

Here, we reveal the details of the structure of the vortex 
sheet by numerically solving Eq. In the following, we 
consider the two-dimensional problem by assuming the 
translation invariance along the z axis {ujzi — 0), reduc- 
ing the original wave function as ^'i(r) = ilJi{x, y)/\/Rz 
with the typical size Rz along the z-direction. It is con- 
venient to measure the length, time and energy scale in 
units of 6ho = \/?i/2mi2(I', and fiiu, respectively, 

with uj ~ {lo\ + aj2)/2. Replacing the wave function 
as — > /h\io, we obtain the dimensionless two- 

dimensional GP equations 



"^12, 



4mi2 



-Q.Lzji^t = lJLiipi. (6) 

Here, we define two-dimensional coupling constants Ui = 
STTaiN.i{mi2 / rrii) / Rz and Ui^^-i = Anai2N3^i/ Rz. Since 
the particle number of each component is conserved, 
the normalization of the wave functions is taken as 
Jdrji/jip = 1. Using the imaginary time propagation 
of the time-dependent version of Eq. ([6]), we calculate 
the equilibrium solutions after sufficient convergence of 
some quantities, such as the total energy, of the system. 

First, for simplicity, we assume mi = m2 = m, 
LUi — UJ2 — Lu, and iVi = N2 = N. The typical example 
of the vortex sheet structure is shown in Fig. [21 where 
we used the parameter values as ui = U2 = u ~ 4000, 
U12/U = 6 = 1.1, and Q. = 0.85a;. The two compo- 
nents undergo phase separation to form several density 
domains of the same component. Goncurrently, as de- 
noted by crosses in Fig. [H^a), the nucleated vortices 
merge to form a winding sheet structure like "serpen- 
tine" instead of forming a periodic lattice. The vortices 
of the -01 component are located at the region of the den- 
sity domains of ip2 component. This can be understood 
from the fact that the condensate density of one compo- 
nent works as a pinning potential for the vortices in the 
other component. By forming vortex sheets, the conden- 
sate achieves remarkable phase separation compared to 
a lattice. 

Figure [2)Ib) shows the radial density and velocity pro- 
file for each component. The domains make a layered 
structure in the radial direction with a certain width. 
The total density profile nx = ni + n2 is nearly smooth 
even in the presence of sharp domain boundaries, being 
given by the Thomas-Fermi profile with 5=1 



UT = \l— — — ( 1 



with the Thomas-Fermi radius 



7?2 



Rtf — 



8u 



1/4 



(7) 



(8) 
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FIG. 2: (Color online) Typical structure of the vortex sheets, 
obtained by the numerical simulation of the coupled Gross- 
Pitaevskii equations. The parameter values are ui = U2 = 
4000, ai2/a = 5 = 1.1, and = 0.85cj (see the text), (a) 
The contour plots of the two-dimensional density profiles of 
the two-component condensates and \tp2'^ ■ The vor- 

tex sheets are made up of chains of singly quantized vortices 
whose positions are marked by x. (b) Solid, dashed, and dot- 
ted curves show the radial density profiles of IV"!!^! l''/'2p, and 
the total density nx = + |V'2|'^, respectively. Circular 

and triangular dots represent the corresponding velocity pro- 
files |vi| = y^l^JiJ^ipH-fSj/Wp with i=l and 2 (in unit of the 
length 6ho and time The velocity of rigid body rotation 

« = f2r is shown for comparison. 

We confirm that for 5 > 1 the total density is always 
approximated by Eq. ([7]), if there is not asymmetry of 
the parameters for intracomponents such as ui ^ u^. On 
the other side, each radial velocity Viir) is flat between 
the vortex sheets, jumping across the sheet to follow the 
velocity of rigid body rotation v = Sir. This behavior is 
quite similar to that in the superfluid helium shown in 

Fig.m 

As b increases, the intercomponent repulsion becomes 
stronger, so that the system favors the smaller region of 
the domain boundary and the larger size of the density 
domain. As a result, the intersheet distance gradually 
increases and the clear serpentine structure disappears 
(Fig. [3]). In the limit of large i5, the equilibrium state 
eventually evolves "rotating droplets," as shown in the 
inset of Fig. [71 where vortex chains do not penetrate into 
the condensate domains and each center of mass of the 
condensate carries the angular momentum. This struc- 
tural change occurs when the equilibrium domain size 
(sheet distance) becomes comparable with the size of the 
condensate (Thomas-Fermi radius), discussed in the next 



subsection. 

We find that there exist different metastable vortex 
sheet configurations with almost the same energies, which 
results from the accidental and inevitable dislocations 
nucleated during the imaginary time propagation from 
different trial initial configurations. As one can see be- 
low, the energy of the vortex-sheet state is mainly de- 
termined by two contributions, kinetic energy of the su- 
perflow and the surface tension, rather than the shape 
of the vortex sheets [111, [l^l . One general feature found 
here is that most vortex sheets are wound. This is in 
contrast to the consideration by Woo et al. 23], who 
expect that the most energetically lowered configuration 
could be the straight vortex sheets (shown in the inset 
of Fig. [7]) instead of the winding sheets. The straight 
sheet configuration is a relatively periodic structure, so 
that this could be obtained near the border of the phase- 
separation regime S ^ 1; for 6 < 1 the stable structure 
is actually periodic vortex lattices [13, HH- For super- 
fiuid He'^ in a rotating container, numerical simulations 
showed that the enertgetically lowest configuration is the 
winding sheets with double spiral pattern, consisting of 
a single piece of sheet with two connections to the wall 

When the intracomponents have asymmetric coupling 
constants ui ^ U2, which corresponds to the system with 
imbalance of the populations, the intracomponent scat- 
tering lengths, or the mass |34{ , the winding vortex sheet 
cannot be obtained, as shown in Fig. |4l One component 
with a smaller coupling constant is confined at the cen- 
ter of the trap, and the other component surrounds like a 
shell [H, [3, [IBl • Since the condensates contain vortices, 
this can be seen as the rotating core-shell structure. In 
Fig. IHa), some density fractions of the ipi component 
are trapped by the vortex cores of the V'2 component. As 
a result, the vortex lattice in the V'2 component possesses 
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(a) (b) (c) 




FIG. 4: (Color online) Typical two-dimensional density pro- 
files and locations of vortices in the vortex states for imbal- 
anced intracomponent parameters; 5 = 1.1, ui = 4000, and 
U2 = 3000 (a), 3500 (b), and 3900 (c). The upper and lower 
panels show IV'iP and |i/'2p, respectively. 



square symmetry due to the antiferromagnetic nature of 
the vortex interactions [21]. With decreasing the imbal- 
ance between ui and U2, the surface of the inner compo- 
nent is distorted to allow the penetration of a chain of 
vortices, as shown in Figs. ID^b) andlDJc). Therefore one 
can expect that the clear vortex sheet can be observable 
for the small parameter imbalance between two compo- 
nents. In experiments, since tuning of both the intra- 
and intercomponent scattering lengths is available by us- 
ing the Feshbach resonance , such imbalance could be 
reduced under the suitable condition. 



C. Estimation of the intersheet distance 

To understand the vortex sheet structure more clearly, 
we give analytical discussions by estimating the inter- 
sheet distance h. According to the numerical result in 
Fig. [2] and the discussion in Sec. [Til we introduce a sim- 
ple planar model for two-component BECs as shown in 
Fig. [5l We assume that there is no asymmetry in the in- 
tracomponent parameters: ui — U2 = u, mi = m2 = m, 
and ui = LU2 — io. Here, the two components form 
straight layered vortex sheets along the y axis, the sheet 
distance b being defined as the length between one sheet 
and its nearest sheet in the same component. This cor- 
responds to the view of the cylindrical planar sheets 
without defects in the polar coordinates instead of the 
Cartesian ones. In equilibrium, the condensate exhibits 
the rigid body rotation V x Vrb — 2riz, following the 
rotation of the reference frame of the system; we take 
Vrb = 2rixy. However, the actual velocity of the conden- 
sate has a staircase like change along the x axis because 
of the vortex sheets. In our model, the velocity of each 
component is assumed to be uniform between the vortex 




FIG. 5: (Color online) Simple model to calculate the inter- 
sheet distance b for the vortex sheet in two-component BECs. 
The rectangular two-component domains align periodically 
and alternately along the x axis with small domain-boundary 
region characterized by the penetration depth \p. The total 
density is uniform all over the space. The condensate veloc- 
ity Vi = Vif has a staircase like increase along the x axis; 
the value jumps 2f26 across the sheet for each component and 
follows the rigid body velocity 2Qx. 

sheets of the same component; the value of Vi increases 
by 2ri6 across every sheet. We neglect the effect of the 
trapping potential, so that the total density nx = ni+n2 
is approximately uniform even in the presence of the do- 
main boundaries with the penetration depth An of the 
order of the healing length = h/^/2rnJl jlj, Isil, Is^ . 
Although the actual condensate density is nonuniform 
because of the trapping potential, this model could be 
valid for condensates with large size in the central region 
of the trapping potential. 

We calculate the energy E of Eq. ^ per volume in 
the range < x < b with the sheet distance b being 
determined so as to minimize E per unit area. The b- 
dependent energy comes from the kinetic energy of the 
condensate flow and the excess energy due to the domain 
boundaries, namely the surface tension a (l6| . The flow 
energy per volume in the rotating frame is 

1 /"^ K - 21^a;)2 1 „ „ 

- / dx } mifii = —mnTll b . (9) 

bJo i 2 48 

Here, the integral was done within the range [6/4,36/4] 
([0, 6/4] and [36/4, 6]) for z = 1 (2) by neglecting the pen- 
etration of the domain. Because of the constant density, 
the interaction energy is independent of the sheet dis- 
tance 6. Then, the 6-dependent energy per unit volume 



6 



is written by 



E 1 ^2:^2 2cr 

— ~ — mnTil b H , 

6 48 b ' 



(10) 



where we included the surface tension from two domain 
boundaries within < x < b. Minimizing this energy 
with respect to 5, one obtains 



6 = 2 



6cr 



1/3 



(11) 



The analytical expression of the surface tension a can be 
found in Refs. (31. , _32]; for mi = m2 = m and ai = 
02 — a, the surface tension in a weak segregated regime 
((5 ~ 1) is given by 



(12) 



and that in a strongly segregated regime ((5 > 1) is given 
by 



4PC 



2V2 0.514 2 /0.055 0.067 



(51/4 



(53/4 J5/4 



(13) 



with the healing length ^ — h/ ^2mfi and the equilibrium 
pressure P — ^\/2gi = ^\/2g2. These two regimes have 
a crossover near 5 ~ 1.34, according to Fig. 3 in Ref. 
[111. In the following, we use Eq. (HH) for (5 < 1.34 and 
Eq. (Uni) for 5 > 1.34 to evaluate Eq. (fTT|) . 

By using the Thomas-Fermi density with (5 = 1 [Eq. 
^] at r = as the value of nx in the denominator of Eq. 
(jlip . one can estimate the approximate sheet distance b 
for the trapped two-component BECs. Figure [6] repre- 
sents the value of 6 as a function of Q. and (5~^. The two 
results using Eqs. and are smoothly connected 
at — 0.75. The sheet distance monotonically de- 
creases (increases) as the rotation (the intercomponent 
repulsion) becomes faster (stronger). We compare our 
numerically estimated values of b taken from Figs. [5] and 
[21 obtaining excellent agreement with Eq. pT|) . 

Besides the sheet distance, we have the other charac- 
teristic length scales in this system: (i) the penetration 
depth Ap of the domain boundary; for the weakly seg- 
regated regime, this is given by ~ S,/VS — 1 (ii) 
the Thomas-Fermi radius Rtf given by Eq. ([8]). By 
comparing these lengths with 6, we can identify three 
distinct structural phases, summarized by the phase dia- 
gram in Fig. [T] The vortex sheet is expected in the region 
b < Rtf- The border b = R^p is gradually increased as 
S goes to infinity, where we can get a critical rotation 
frequency flc = [16'K/{y/mi + 167r)]i/^ above which the 
vortex sheet is observable for arbitary 5. In the region 
b > i?TF , the vortex sheet cannot penetrate into the con- 
densate domain, the condensates forming droplets with a 
center-of-mass rotation [see the inset (a) of Fig. [7]. Since 
most of the vortices are located in the low-density region 
as "ghosts" [S^, they contribute little to the energy cost. 
Thus, the rotating droplet can exist even for very small 

n. 
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FIG. 6: (Color online) The sheet distance b (in unit of 6ho) 
as a function of Q, (a) and (b) for u = 4000. In (b), the 
two results using Eqs. (|12p and (I13|l are smoothly connected 
at = 0.75. Also, the numerically estimated values of h for 
Q. = 0.85 are plotted for 5 = 1.1 (Fig. 0, 1.5, 2.0, and 3.0 
(Fig. [3)l; the values were taken by averaging the intersheet 
distance from the profile of vortex positions. 



In the region b < Rtf, the equilibrium structure is 
generally the serpentine vortex sheets. However, when 
b becomes comparable to the penetration depth Ap, the 
periodic stripe structure of vortex sheets is formed [Fig. 
Wl,c)] instead of the serpentine one. We determine this 
border at 6 = 4Ap, where the density plateau for each 
component is vanished in our model of Fig. O In this 
weakly segregated regime, the influence of the periodicity 
of vortex states for S < 1 survives. As b increases with 
(5, the domain size increases, which leads to the unfavor- 
able tight confining of the relevant number of vortices 
into the decreasing low-density regions. As a result, the 
stripe configuration is disrupted by linking edges of the 
neighboring sheets; the inset (b) in Fig. [7] is obtained 
by the imaginary time propagation and the gradual in- 
crease in (5 from the state (c). Therefore a tendency to 
confine vortices as loosely as possible and the inevitable 
dislocations give rise to the serpentine vortex sheets such 
as Figs. [5] and [31 although the total energy of the qua- 
sistripe state (c) is lower than those configurations. This 
mechanism is involved when the vortex sheets is created 
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Vortex stripe sheet 




FIG. 7: (Color online) Q-S~^ phase diagram of the equilibrium 
vortex states. The borders are determined by comparing the 
characteristic length scales (see the text). The insets show 
the typical equilibrium structures of the density Itpif and 
the vortex positions for the parameters indicated by the star 
marks; (a) {Q,5) = (0.4,3.0), (b) (^1,5) = (0.85,1.5), (a) 
(f2, 5) — (0.85, 1.02). In (b), the region indicated by a circle 
represents the linked edges of two vortex sheets. 



If the imbalance of the intracomponent parameter is suf- 
ficiently small, the two-component BECs form a winding 
vortex sheet structure; otherwise, the condensates form 
a rotating core-shell structure. We developed a simple 
model to estimate the vortex sheet structure, following 
the argument of superfluid helium by Landau and Lif- 
shitz. As the intercomponent coupling becomes larger, 
the sheet distance is also increased, approaching the ra- 
dius of the condensate density. With this increase, the 
sheet structure changes from the stripe pattern to the 
winding one like "serpentine." In the slowly rotating case, 
the sheet distance exceeds the condensate radius so that 
the condensate forms rotating droplets, where the angu- 
lar momentum was carried by their center-of-mass rota- 
tion. These results are summarized in the phase diagram 
of Fig. [71 obtained by the comparison of the character- 
istic length scales, the sheet distance b, the penetration 
depth Ap, and the Thomas- Fermi radius Rtf- 

Two-component BECs with tunable intercomponent 
interactions were obtained experimentally p^ . [llj . By 
rotating them, it is possible to explore the interesting 
regime of vortex phases in multicomponent BECs. 
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